Magnetic breakdown in ortho-II high-temperature cuprates 
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The recent confirmation of closed Fermi surface in the high temperature cuprates via de Haas- van 
Alphen and Shubnikov-de haas oscillations calls for a theoretical investigation of the origin of such 
oscillations. We study de Haas-van Alphen oscillations and magnetic breakdown in the context 
of ortho-II high temperature cuprates to understand the origin of the multiple frequencies. We 
find that the magnetic breakdown is highly sensitive not only to the ortho-II potential, but also to 
the Fermi surface topology, and is thus useful to distinguish between various theoretical proposals 
related to quantum oscillations. 
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I. INTRODUCTION 



A remarkable phenomenon associated with a Fermi 
surface in a metal is the observation of periodic oscil- 
lations in thermodynamic and transport quantities as a 
function of the inverse of the magnetic field. This quan- 
tum effect is due to a quantization of Landau levels, and 
its periodicity is proportional to the area enclosed by the 
Fermi surface in two dimensions (or the area enclosed 
by extremal orbits in three dimensions). One of the un- 
solved puzzles in the high temperature (high Tc) cuprates 
is a truncated Fermi surface, dubbed the Fermi arc, 
detected by angle resolved photoemission spectroscopy 
(ARPES)ii"— This cannot be understood within a con- 
ventional theory - for instance, if the truncation occurred 
due to a broken translational symmetry, this would lead 
to Brillouin zone folding, which has not been observed so 
far in ARPESi^ This unconventional metallic state with 
a Fermi arc has been called the pseudogap phase. Up 
until 2007, it seemed that quantum oscillations failed to 
occur in the pseudogap state of high Tc cuprates. 

Since the first striking report of the high Tc 
cuprates by Doiron-Leyraud et al^ in 2007, a series 
of experiments have reported quantum oscillations in 
YBa2Cu306.5i(YBC06.5i) and YBa2Cu408, high purity 
high Tc materialsii^— The main frequency of oscillation 
of 540T in YBCOg.si corresponds to 2 % of the area of 
the Brillouin zone. This area deviates dramatically from 
the nominal doping of 10% expected from the Luttinger 
sum rule. In contrast, overdoped Tl2Ba2Cu06+2; shows 
a single hole Fermi surface^^ which correspond to 1 + p 
holes per Cu satisfying the sum rule and further con- 
firmed by the Hall coefficient RjyJ^ A recent report dis- 
played a similar contrast between overdoped samples and 
underdoped ones in the electron-doped Nd2-xCe2;Cu04 
(NCCO) where a large electron-like pocket is reported 
at X = 0.17, whereas small hole-like pockets are re- 
ported at a; = 0.15 and x = 0.16^^ Shortly after the 
experimental discoveries, several theoretical proposals 
have been made J^"— Most of them share the underlying 
idea that the quantum oscillation arise due to a broken 



symmetryji^Ti^^ Some have further offered a missing link 
between the Fermi arc and quantum oscillation from a 
closed Fermi surfacej^i^ 

However, the precise nature of the broken symmetry 
is still under debatei^^r^l Can we identify the order by 
analyzing the current experimental data of quantum os- 
cillations? One salient aspect of the available experimen- 
tal data is a discrepancy between experimental reports in 
YBCOg.si- While they all agree that the main frequency 
of oscillation is about 540 T, the satellite frequencies are 
different. LeBoeuf et alX confirmed the first observed 
main frequency of 540 Tj^ In addition, they found satel- 
lite frequencies of 450 T, 630 T (suggested to be due to a 
combination of warping of the Fermi surface and bilayer 
splitting), and 1130 T (suggested to be a second harmonic 
of the 540T main frequency) . On the other hand, whereas 
Sebastian et alr^ found the main frequency of 540 T, they 
also observed an additional satellite frequency of 1650 T. 
Recently, Riggs et al. reported similar oscillation fre- 
quencies in specific heat measurements as those found 
by Sebastian et alM We note that multiple frequencies 
have so far been found only in ortho-II material and its 
understanding requires a study of magnetic breakdown. 

The ortho-II potential present in YBCOe.si is not a 
strong potential, yet it is visible in some experimental 
probes. For example, while the effect of the ortho-II po- 
tential on ARPES is negligible^, the potential has been 
detected in Raman spectroscopy^^. One may naively ex- 
pect that such a potential can be screened as pointed out 
in Ref. [2^ However, first principle studies show that the 
oxygen ordering leads to order in the Ba positioii^, which 
leads to a stronger ortho-II potential than one would 
otherwise expect. A Fermi-surface-induced lattice mod- 
ulation was observed only in oxygen empty CuO chains 
via diffusive x-ray scattering measurements^ indicating 
the important role of oxygen ordering in electronic struc- 
tures. It was recently suggested that the ortho-II po- 
tential can be used to differentiate between two (tt, tt) 
orders, antiferromagnetism and d-density wave ordersi^ 
In this scenario, the ortho-II potential plays an impor- 
tant role, but its relative weakness implies that magnetic 
breakdown effects must be taken into account explicitly. 
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For instance, if the magnetic breakdown occurs below the 
field at which oscillations begin to appear (30 T), then 
the ortho-II potential can be neglected altogether in the 
analysis of quantum oscillations. On the other hand, if 
the breakdown occurs within the currently accessible field 
range of 30-60T, then two different sets of Fermi surfaces 
- corresponding to the Fermi surface before and after the 
breakdown - may participate in the quantum oscillation. 
In this case, the dominant oscillation frequencies would 
depend sensitively on the field range at which a particular 
experiment is carried out. 

In this paper, we ofi^er a general discussion on quantum 
oscillations in Sec. |ll]and define a way to take magnetic 
breakdown into account for non-trivial Fermi surface ge- 
ometries. In Sec. mil we focus on the AF-|-ortho-II sys- 
tem for its relevance to the ortho-II materials studied 
experimentally. In Sec. IIVI we study incommensurate 
orders, and show the Fermi surface of spin spiral order 
in the presence of ortho-II potential. We determine os- 
cillation frequencies and magnetic breakdown fields, and 
compare the result with those obtained for the commen- 
surate order. We summarize our results and discuss the 
implications of our findings in Sec. |Vl 



II. MAGNETIC BREAKDOWN AND FERMI 
SURFACE CURVATURE 



Quantum oscillations arc a quantum mechanical effect 
observed in metals with closed Fermi surfaces at low tem- 
peratures in the presence of a magnetic field. As the 
intensity of the applied field increases, there are oscilla- 
tions in the physical properties of the metal, including the 
magnetic moment (named de Haas- van Alphen, dHvA), 
resistivity (Shubnikov-de Haas, SdH), specific heat, and 
sound attenuation. When plotted versus the inverse field, 
these quantities display a remarkable regularity in their 
oscillations. 

Quantum oscillations can be understood via the semi- 
classical quantization of quasiparticle energies in an ap- 
plied magnetic field. They are of great significance for 
their relation with Fermi surface geometry: the frequency 
of oscillations (in Tcsla) is related to the area enclosed by 
a Fermi pocket (in 2 dimensions) or the area enclosed by 
an extremal orbit of the Fermi surface (in 3 dimensions) 
by the following equation: 
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where Aj- is the area enclosed by the Fermi pocket and 
F is the frequency of oscillations. If there are multiple 
pockets (or multiple extremal orbits), that will result in 
multiple frequencies that can be determined by Fourier 
transform. 

The different frequencies are constrained by the Lut- 
tinger sum rule which relates the density of carriers to 



the area enclosed by the Fermi pocket through: 
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where a, b are the lattice constants and the factor of 2 
comes from spin degree of freedom. Therefore, we can 
relate the carrier density to the frequency of oscillation 
found through: 



(3) 



where $o is the unit quantum fiux. One must be careful 
when working with a broken translational symmetry to 
work in the reduced Brillouin zone (RBZ). In the cases 
that we will study, this will be important. The Luttingcr 
sum rule specifies that the total doping p is equal to the 
sum of carrier density for each pocket, that is 
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where h.l.p mean hole-like pocket and e.l.p means 
electron-like pocket. We can easily convert this equation 
using the frequency of each pocket: 



i—h.l.p j— e.l.p 



(5) 



This equation establishes a constraint on the frequencies 
of the system. 

Magnetic breakdown is an effect that occurs at large 
magnetic fields for systems with multiple Fermi pock- 
ets that almost intersect one-another. As first discussed 
by Cohen and Faliko\«^, magnetic breakdown happens 
when the transition amplitude for the quasiparticle to 
tunnel from one band to the other become sizeable as 
the magnetic field increases. A semiclassical treatment 
shows that magnetic breakdown occurs approximately 
when the condition TilJc^f > E'^ is satisfied^. Here, ep 
is the Fermi energy, and Eg is the energy splitting of two 
bands generated by a small perturbation such as a spin 
density wave, a charge density wave order or some lattice 
potential. In terms of the breakdown field (-B*), this can 
be written as 
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A more general treatment^i^ yields a useful formula 
involving only the local geometry of the Fermi surface to 
find the breakdown field which, in the case of almost free 
electrons reduces to 



B* = 
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where 5k is the k-space distance between the two bands 
at the chemical potential in units of 1/a. This equation 
can also be written as 
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where £b = ^/<^o/B*. The constant of V2 in Refs^i^s 
is obtained using a circular orbit as the local geometry 
of the Fermi surface. However in general, this constant 
will depend on the details of the geometry of the Fermi 
surface at the band separation. We therefore define a 
constant K as 



K. 



(9) 



Below, we will show how K depends on the FS curvature. 

It is important to note that magnetic breakdown does 
not occur at a sharply-defined field, but is a crossover 
effect. It is usually defined to be the field at which the 
probability for tunnelling is equal to exp{~l). In our 
numerical simulations below, we will define the magnetic 
breakdown field to be the field at which the intensity 
of oscillations arising from the pockets in the original 
Fermi surface is equal to the intensity of the pockets af- 
ter reconstruction of the Fermi surface. Therefore, one 
should note that the intensity of oscillations from the re- 
constructed Fermi surface is still visible after the break- 
down field as defined in this paper. 

To study the effect of curvature on magnetic break- 
down, we compute the value of K for various systems 
that have different curvature. Our numerical simulations 
consist of solving a lattice model in an applied magnetic 
field. The mean-field Hamiltonian is given by 

ija i 

+AoJ2{-iy (cl^c^t - 4i^a) ■ (10) 



Here, c]^ creates an electron with spin a at site i, is 
the chemical potential, A is the ortho-II potential, which 
is staggered in the x-dircction, and Aq is the order pa- 
rameter of the broken symmetry, like an antifcrromag- 
net. The amplitudes tij are real, and we keep nearest 
neighbor (t), next-nearest neighbor (i'), and third near- 
est neighbor (t") hoppings. In addition, the orbital effect 
of magnetic field has been introduced in the standard 
manner through the Peierls substitution ty — > ty-e'^*^, 
where Aij = ^ J^' di ■ A, and A is the vector potential. 

The three specific systems used to study the Fermi 
surface curvature dependence of K are (a): A finite and 
Ao = , (b): A = and Aq finite and (c): both A and 
Aq finite. The Fermi surface of these three systems arc 
displayed in Fig. [TJ 

We obtain the quasiparticle spectrum of the total 
Hamiltonian by exact diagonalization of Bloch wave func- 
tions on a magnetic unit cell of size x The in- 
clusion of an orbital magnetic field implies that, for any 
given fixed gauge, the vector potential appears to break 
translational invariance. Translational invariance is re- 
stored to the system by noting that translations by 
and Ly must be accompanied by appropriate gauge trans- 
formations. This can only be done in a consistent fashion 
provided that the magnetic fiux piercing the magnetic 





FIG. 1: Fermi Surfaces for three different systems studied to 
explore magnetic breakdown systematically. The parameters 
for the electronic structure are chosen to generate different 
Fermi surface curvature to study its effect on magnetic break- 
down. Note that the figures show only one quadrant of the 
Brillouin zone. 



unit cell is quantized, L^LyB — n^Q^ Here n is an in- 
teger, and $0 = hc/e is the quantum of flux. In our 
simulations, we fix the height of the magnetic unit cell 
Ly = 2, and also fix a single fiux quantum per magnetic 
unit cell. We then sweep the magnetic field by varying 
the width of the magnetic unit cell, 1/B = 2Lx/^o- 

In each system, in order to compute K, we vary the 
different parameters in order to obtain different values 
of Sk- For each value of 5k, we then compute the en- 
ergy, the density of states (DOS) and magnetization of 
the system as a function of magnetic field. Fig. [2] shows 
an example of magnetization as a function of the inverse 
field. We then compute the breakdown field by use of 
a "field-dependent Fourier transform" , which consists of 
computing the Fourier spectrum over a narrow window 
of fields, and shifting the position of these windows to 
get a running field dependence for the Fourier spectrum. 
We require the window of fields to be broad enough to 
contain a few oscillations of all relevant frequencies, but 
narrow enough to detect the magnetic breakdown. We 
then compare the intensities of the frequencies from the 
band present before the breakdown to those after the 
breakdown and determine the breakdown field to be the 
one at which the intensities are equal. Fig. [3] shows this 
comparison for the example of Fig. [21 From this graph, 
we determine B* to be approximately 25T. We compute 
Ib from B* , and obtain 5^ by determining the k-space 
separation between different pockets at the chemical po- 
tential; K is then simply the product of Ib and 5k- 

l(a)[ the parame- 



In system (a). 



shown in Fig. 
ters used, written in the (t, i', t", /i. A, Aq) basis, were 
(0.3,-0.09,0.012,-0.25, A, 0). We tuned 4 by vary- 
ing A over the values A = 0.022, 0.025, 0.03, 0.035, 
and 0.05. The resulting value of K is shown as the 
dashed line in Fig. |4]and is approximately 1.8. For sys- 
tem (b), shown in Fig. l(b)[ we chose the parameters 
(0.3, —0.06, 0, —0.5, 0, Ao) and we tuned 5k by varying the 
value of Ao over the values Ao = 0.03, 0.035, and 0.04. 
The resulting value of K is shown as the dotted line in 
Fig. [4] and is appro ximately 3.5. Finally for system (c), 
shown in Fig. 1(c) the parameters used, written in the 
(t,i',i",^,A, A^^basis, were (0.3, -0.09, 0.012, -0.27, 
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FIG. 2: de Haas-van Alphen effect for the AF + 
ortho-II simulation. The parameters used for this 
case, written in the {t,t' ,t" , fi, X, Ao) basis, were 
(0.3, -0.09, 0.012, -0.27, 0.08, 0.07) 
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FIG. 3: Field-dependent Fourier intensities of the /3 and 7 
pockets summed together (solid line), compared to the inten- 
sity of the P' pocket (dashed line), for the data in Fig. [2] 
From this graph we estimate the breakdown field B* 25 T. 



0.08, 0.07), (0.3, -0.095, 0.008, -0.271, 0.07, 0.08), (0.3, 
-0.09, 0.012, -0.27, 0.1, 0.07) and (0.3, -0.09, 0.01, 
-0.27, 0.05, 0.21), giving four different values of 4- The 
resulting value of K is shown as the solid line in Fig. |4] 
and is approximately 2.2. Our results confirm that B* 
strongly depends on the Fermi surface curvature. Note 
that when K is doubled, B* is quadrupled. 



III. QUANTUM OSCILLATIONS IN ORTHO-II 
+ AF SYSTEM 

From the three systems studied in the previous sec- 
tion, case (c) is the one that is relevant to the ortho-II 
YBCO materials. The different Fermi pockets for this 
system are shown and labelled in Fig. [SJ The K param- 
eter for this system was found to be around 2.2, which 
differs from for the circular Fermi surface^^^. Since 
B* oc K^, this difference leads to a breakdown field that 
is two times larger than what one would get from a naive 
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FIG. 4: K parameter as a function of Sk- We first note the K 
is indeed roughly constant for a given curvature. The dashed 
line is only ortho-II (case (a)) and the dotted line is only 
AF (case (b)). The solid line represents the AF -|- ortho-II 
system (case (c)). Note that K is biggest for case B where 
the curvature of Fermi surface is sharp, so a weaker magnetic 
breakdown field is expected. 




FIG. 5: The different Fermi pockets of the AF-|-ortho-II 
and DDW-|-ortho-II potentials. In a single layer system, the 
AF-|-ortho-II surface has a, /3 and 7 pockets whereas the 
DDW+ortho-II surface has a and /?' pockets. 

magnetic breakdown analysis. For the parameters used 
in last section with Aq = 0.07, we expect oscillations at 
Fa ^ 524 T (electron pocket), Ffj ~ 1630 T (hole pocket) 
and ^ 328 T (hole pocket). The smallest band split- 
ting occurs between the /3 and 7 bands. Hence, at large 
magnetic fields, magnetic breakdown will lead to a com- 
bination of the /3 and 7 pockets into /?' pockets, with a 
frequency that is approximately equal to the average of 
the 13 and 7 frequencies, F'^ - ^^^°+^^'^ = 979 T. From 
the magnetization shown in Fig. [5J a Fourier transform 
yields the spectrum shown in Fig. [S] where all the ex- 
pected frequencies can be seen. 

To compute the breakdown field in this case, we can 
make use of the /C-parameter computed in the previous 
section. The only thing we need in order to find B* is 6k 
which can be computed analytically in this case from the 
dispersion of the /3 and 7 along the ky = Tr/2 line. These 
dispersions are given by 

= -2t"(cos 2k^ - 1) - y/At^ cos2 k^ + {A± A)2. (11) 

Note that Sk is equal to the difference between the val- 
ues of k^ at which i?^(fc^) = /i. t, t", A, and A 
have a direct effect on Sk, while t' docs not. Increas- 
ing t and t" lowers Sk- It is straightforward to see 
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FIG. 6: Fourier transform of the dHvA data shown in Fig. [2l 
This spectrum shows the following peaks; a, /3, 7, /?', /?' + 
7,/3' + /3. The parameters are {t,t' ,t" , n, X, Aq) — 
(0.3,-0.09,0.012,-0.27,0.08,0.07). Here we used a 
Lorentzian broadening of 10^'^ for Landau levels, and the field 
range for the transform was 20 - 100 T 



broaden the Landau levels from delta functions in energy 
to Lorentzians of width F w l/r. For instance, in the 
previous section, we chose F = 1 meV— , but F may differ 
between samples. Here we study how varying the value 
of F affects the Fourier spectrum. Table U shows the 
Fourier intensities, computed in the range of 20 — 70 T, 
of all the pockets for different broadening. Note that the 
oscillation intensity is more sensitive to disorder for some 
pockets than for others; hence, broadening has a direct 
effect on B*, which is also shown in the table. 



TABLE L Fourier spectrum as a function of broadening. 
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the effect if one takes either t or t" to 0. For exam- 
ple, when t" = 0, Sk = arccos(^M^ - (A + A)2/2t) - 
arccos(-\/A'2 — (A — X)'^/2t). If the band is normalized 
and the effective band mass is bigger (effectively lowering 
t), it will push B* higher. On the other hand, increas- 
ing A and A increases Sk, as we expect. While t' does 
not have a direct effect, it affects Sk via shifting of the 
chemical potential to adjust for a constant density. 

While we obtain B* ^ 25T for the set of parame- 
ters given above using K = 2.2, we emphasize that Sk 
is very sensitive to the choice of parameters discussed 
above. For example, B* - lOOT for A = 0.1 and A = 0.1 
{6k ~ 0.13). Eq. (|lip is only true for antiferromagnetism. 
Note that an ortho-II system with DDW ordering has a 
similar zone folding to the AF+Ortho-II case considered 
here. However, as discussed in Ref.[2^ in the DDW case, 
symmetries prevent the splitting of the /3 and 7 bands. 
Hence, for DDW order in the presence of ortho-II poten- 
tial (DDW-|-ortho-II order) in monolayer systems, there 
are only a and /?' pockets and no magnetic breakdown 
occurs within the accessible field range. On the other 
hand, in bilayer systems such as ortho-II YBCO, the (3 
and 7 bands do split if there is a current circulating be- 
tween layers as shown in Ref. [22I . Thus, simulations of 
AF-f-ortho-H on a monolayer give us direct information 
regarding magnetic breakdown of DDW+ortho-II order 
on a bilayer, without the additional computing resources 
necessary to simulate the bilayer system. 



A. Disorder effect and anisotropic scattering 

The different ortho-II materials studied so far possess 
a rather large mean-free path of order / « 160 Ai^. Using 
an average Fermi velocity of 8.4 x 10** m/si^, we get a 
scattering rate l/r « i;^?/? ~ 3 meV. In order to simulate 
the effects of disorder that exist in real materials, we 



Note that for F = 2 meV, the oscillation intensity from 
the /3 pocket, is still visible, but that for F = 3 meV, 
I^ is gone. Hence a small difference in F can make a qual- 
itative difference in the Fourier spectrum. This is under- 
standable because the broadening of the Landau levels 
will affect the larger frequencies more than the smaller 
ones. This is why I^ goes from being the smallest at 
F 0.6 meV to being the largest at F = 3 meV (see 
Table H]). The 7 pocket, having the smallest frequency, is 
also the least affected by broadening. 

Another important aspect of Fourier spectrum pro- 
duced from our analysis is the difference between the 
relative intensities of our multiple peaks and those ob- 
served in experiment. For example, in our case, the /?' 
pocket is the dominant peak unless the disorder broad- 
ening is greater than 3 meV, while its measured ampli- 
tude is small, and not confirmed by all of the experi- 
mental groups. One way to reconcile the difference is 
to introduce an anisotropic scattering on the Fermi sur- 
face. If the quasiparticle scattering rate is momentum- 
dependent, then it is possible to obtain the electron 
pocket, a, to be the dominant frequency. For example, if 
quasiparticles near k = (7r/2,7r/2) (where the hole pock- 
ets are) have a larger scattering rate than the quasipar- 
ticles near k = (tt, 0) and k = (0, tt) (where the electron 
pockets are), then the a pocket could get the largest 
intensity. However, this contrasts with the anisotropic 
scattering observed by ARPES at high temperatures in 
the absence of a magnetic field. It was shown that the 
scattering rate deduced from the imaginary part of the 
self energy is lower along the node compared to the antin- 
ode directiori^^. However, as we discussed in the intro- 
duction, the ARPES have shown the Fermi arc at high 
temperatures (above Tc) without the magnetic field, and 
it is possible that the low temperature and high field state 
may not be smoothly connected to a high temperature 
and zero field pseudogap phase. 
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FIG. 7: Fermi surface for spiral spin density wave order at 
wave vector Q = (77r/8,7r) in an ortho-II potential. In ad- 
dition, the spiral state, described by Eq. (|12[1. gives rise to 
a second set of Fermi surfaces, not shown here, obtained by 
setting Q — > — Q and flipping up and down spins. 



FIG. 8: Pockets arising from the Fermi surface in Fig.[7l Os- 
cillations at the electron hole frequency Fa ~ 525 T should 
be seen at all experimentally-relevant values of field. On 
the other hand, the hole pockets undergo breakdown at two 
different values of field. For B < B^, oscillations occur 
at Ffi ~ 1650 T and Fj ~ 278 T. At intermediate fields, 
Bl < B < B2, hole pocket oscillations are dominated by the 
frequency Fjs' ~ 965 T. For B > B2, oscillations occur at 
Fpn ~ 1505 T and F^, ~ 425 T. Based on the value of Sk 
responsible for the two breakdowns, we expect B2 / Bl ~ 5. 



IV. INCOMMENSURATE ORDERS 

In both La2_2;Sr2;Cu04 and YBa2Cu306-f-2;, neutron 
scattering experiments have reported strong inelastic sig- 
nals at incommensurate wave vectors, (tt — 5, 7r) i^^i'^° 
Later, it was found that the quasi-elastic incommensu- 
rate peak intensity increases linearly with magnetic field 
in YBC06.45'^ It was also proposed that an incommen- 
surate spiral order is induced by a magnetic field lead- 
ing to a Fermi surface reconstruction responsible for the 
quantum oscillations It is worthwhile to investigate the 
effect of the ortho-II potential for the incommensurate 
spiral order case. While our analysis is specific to the spi- 
ral spin density wave, it can be generalized to a coUinear 
incommensurate case if the higher order gaps are smaller 
than the gap generated by the ortho-II potential42ri^ In 
other words, if we keep only the first order gap generated 
by an incommensurate collinear order our analysis leads 
to the same results as for the spiral case. An incommen- 
surate spiral order breaks translational symmetry just as 
collinear incommensurate order does, but it mixes |fc,t) 
with \k -\- Q, I) and \k, |) with |fc — Q, f). The mean field 
interaction term can be written as Eq. (|12p . 

-ffspiral = Asp ^ cj^^Ck+QJ. + h.C. (12) 
k 

The Fermi surface for the choice of parameters t = 0.3, 
t' = -0.09, t" ^ 0.012, n = -0.27, A = 0.05, and A^p = 
0.08, and for the wave vector Q = (77r/8,7r) is shown in 
Fig.H 

From Fig. [71 there are now two Fermi surface splittings 
between /3 and the 7 pockets which may be relevant to 
magnetic breakdown. We can determine the two different 
band separation at the chemical potential (6k) that will 
then give us an estimate for the two breakdown field 
and i?2 (see fig. [51 They are respectively 6ki = 0.07 and 
6k2 = 0.19. Assuming that K is of order 2 (a reasonable 
estimate since the curvatures at the two separations are 



very similar to the commensurate case), then we get the 
following values for the two breakdown fields: = 35T 
and B^ = 251T. 

V. SUMMARY AND DISCUSSION 

ARPES and dHvA/SdH oscillations arc two of the 
most powerful tools to measure the Fermi surface of a 
metal. The two methods give complimentary informa- 
tion about a system. Quantum oscillation experiments 
probe the bulk of a system, and give a measurement of 
the size of various Fermi pockets in the system. However, 
they reveal neither the location of the pockets, nor the 
number of different pockets of equal area. In addition, 
quantum oscillation experiments must be carried out in 
a magnetic field. In contrast, ARPES, which can give 
details of the locations and shapes of different pockets 
in the Brillouin zone, is a surface probe which must be 
carried out in the absence of a magnetic field. There- 
fore, a combination of the two techniques should offer 
a relatively complete picture of the Fermi surface of a 
material. 

In the high Tc cupratcs, there is a qualitative discrep- 
ancy between the results provided by the two measure- 
ments: ARPES measurements indicate Fermi arcs, while 
quantum oscillations results reveal the existence of closed 
Fermi surfaces. Several theoretica l^^i^" and experimental 
proposals support an electron pocket based on the nega- 
tive Hall coefficient measured at low temperatures^ and 
more recently, high-field thermoelectric measurements 
performed on a series of high Tc materials also indi- 
cate that this 540T frequency comes from an electron- 
like pockeli^. However, there is no visible spectral weight 
such as an electron pocket near (tt, 0) in ARPES. There is 
also a discrepancy in the presence of satellite frequencies, 
such as 1650T, among the different quantum oscillation 
experiments and related proposals. 

Another important experimental observation is the 
absence of a visible Zeeman effect on the phase of 



7 



oscillation4i Sebastian et al proposed that the order re- 
sponsible for quantum oscillations is in the spin triplet 
channel. While the DDW discussed above is in the sin- 
glet channel, it is interesting to note that nematicity al- 
lows a coupling between spin triplet DDW and AF order 
in the presence of a magnetic field.— The nematicity, 
a broken x-y symmetry is reported in YBCO by neu- 
tron scattering!^ Such a coupling leads to a magnetic 
field-induced ordering which is a combination of AF and 
triplet DDW, where the dominant order is the one with 
the smaller gap in the spectruni^ 

In this paper, motivated by a series of quantum oscil- 
lation experiments performed on ortho-II YBCOe.si at 
low temperatures and high magnetic fields, we investi- 
gated the magnetic breakdown effect and the Fermi sur- 
face topology in ortho-II high Tc cupratcs and showed 
how to reconcile the discrepancy in the observed satellite 
frequencies. Our results also pose constraints in a theory 
of quantum oscillations in high cuprates in general. 
We provided a general criterion for the magnetic break- 
down, showing that it depends, in addition to a band 
separation and the Fermi velocity derived from the semi- 
classical study, on the Fermi surface curvature. 

Applying our results to (tt, tt) orders such as AF in the 
presence of ortho-II potential, we found that the mag- 
netic breakdown field is highly sensitive to the ortho-II 
potential, the electronic dispersion, and the order pa- 
rameter strength. As expected, the magnetic breakdown 
field increases when the potential increases, the AF or- 
der get stronger, or the electronic dispersion gets fiattcr. 
Smoother Fermi surface curvature, modified by the po- 
tential, also increases the breakdown field. We found 
that the magnetic breakdown field (defined as the field 
strength at which the intensities of the original and mod- 
ified frequencies match) can change from 25 T to 100 T. 
For example, a factor of 4 increase can be achieved if the 
ortho-II potential (A) increases from 0.08 to 0.1 and the 
AF order strength (A) from 0.07 to 0.1. DDW in bilay- 
ers gives similar results, if there is a circulating current 
between the layers. 

While the shape of the Fermi surface depends on an 
exact broken symmetry (as shown in Fig. 2 for the AF -f 
ortho-II case and in Fig. 10 for the incommensurate spi- 



ral spin density wave), the presence of multiple frequen- 
cies strongly supports a broken translational symmetry. 
Among several frequencies observed in experiments, the 
presence of the 1650T frequency can be explained by a 
/3-pocket that occurs due to a near (tt, 7r)-folding and a 
(tt, 0)-folding (from ortho-II potential) of the Fermi sur- 
face and that such pocket stop yielding oscillations be- 
yond B*. 

In our analysis, the intensity of each frequency depends 
on the disorder strength. A disorder broadening F big- 
ger than 3 meV (which is consistent with experimental 
values) makes the (3 pocket, corresponding to the 1650T 
frequency, disappear. This is simply because a larger 
broadening makes faster oscillations (larger frequency) 
less visible. The disorder strength also affects the inten- 
sity of the other frequencies. As demonstrated in Table I, 
the 7 pocket has the highest intensity for a 3 meV broad- 
ening. One could argue that the main frequency is a hole 
pocket (like the 7 pocket in our model)^. However, the 
difficulty with that argument is that a main frequency of 
660 T was observed in YBa2Cu408, an ortho-II free ma- 
terial. As pointed out in Rcf. [l^, the electron pocket is 
insensitive to the ortho-II potential while the hole pocket 
is strongly modified. Therefore if a similar small fre- 
quency is a universal phenomenon in high Tc cuprates, 
the main frequency should be one that is not modified 
by the ortho-II potential, such as the a-pocket in this 
model. 

In addition, we extended our study to systems with 
incommensurate orders, which undergo a series of break- 
downs with distinctly different Fermi surface shapes. To 
make progress in understanding the topology of the Fermi 
surface and associated order, a systematic study of the 
magnetic field angle dependence on the quantum oscil- 
lations and further ARPES studies in both ortho-II and 
ortho-II-free materials are necessary. 
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